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CALCULATION OF THREE-DIMENSIONAL UNSTEADY TRANSONIC FLOWS 

PAST HELICOPTER BLADES 

J. J. Chattot* 

Ames Research Center 
and 

Aeromechanics Laboratory 
AVRADCOM Research and Technology Laboratories 


SUMMARY 


A finite difference code for predicting the high-speed flow over the 
advancing helicopter rotor is presented. The code solves the low-frequency, 
transonic small disturbance equation and is suitable for modeling the effects 
of advancing blade unsteadiness on blades of nearly arbitrary planform. The 
method employs a quasi-conservative mixed differencing scheme and solves the 
resulting difference equations by an alternating direction scheme. Computed 
results show good agreement with experimental blade pressure data and illus- 
trate some of the effects of varying the rotor planform. The flow unsteadiness 
is shown to be an indispensible part of a transonic solution. It is also shown 
that, close to the tip at high advance ratio, cross-flow effects can signifi- 
cantly affect the solution. 


INTRODUCTION 


Air flow past a helicopter rotor blade exhibits many very complex features 
such as three-dimensional unsteady effects, shock-wave motions, vortex inter- 
actions, and stall. A complete numerical simulation cannot even be attempted 
yet, but it is possible with the present-day computers and numerical methods 
to model some of these features and acquire a better understanding of some of 
the mechanisms involved. 

The model used in this study is a perfect fluid model that is further 
simplified by the small-disturbance approximation. Weak, almost normal shock 
waves are accounted for by retaining the leading nonlinear term in the stream- 
wise direction. This model is useful for simulating the subsonic and transonic 
flow past the advancing blade. Under these conditions the incidence is usually 
small, and the results presented correspond to nonlifting blades. A proper 
wake representation is required to extend this simulation to lifting configur- 
ations. Prediction of the complicated rotor vortex structure is not within 
the scope of the present work. 
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EQUATION AND BOUNDARY CONDITIONS 


The mathematical model used in this report is the three-dimensional 
unsteady (low-frequency) small-disturbance transonic equation as derived by 
M. P. Isom (ref. 1, p. 20). This equation is derived in a blade-attached 
Cartesian coordinate system under the usual assumptions: 

1 - M 2 (l + y) 2 = 0(6 2/3 ) 


where 



R 


c 

Q, 


V 


e = 0(6) 

tip Mach number due to the blade rotation 

advance ratio 

blade thickness 

inverse of the aspect ratio 

blade radius 
chord of reference 
rotational velocity 
sound speed 

forward velocity of the rotor 
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In condensed notation the equation can be written: 


where 


3 2 <t> 
3t 3x 


-lr[ B 



+ 


C 


3 2 j> 
3x 3y 


+ D 


3 2 <t> 

3y 2 


+ E 


3 2 ji 

3z 2 


( 1 ) 


a “ 2mZ (y + ^ cos 

i-j _ 1 - H 2 (y + y cos t ) 2 
B ' {*/> 

B' ■ ^ ^ ^ M 2 (y + v cos t) 

C = 2M 2 y sin t(y + y cos t) 



E = 1 


t, x, y, and z are the dimensionless dependent variables normalized by ft, 
c, R, and 6” 1 ' 3 c, respectively, and y is the ratio of the specific heats. 

At each time step NS in THREED, the coefficients are computed and stored in 
one- dimensional arrays 

A(J) , B(J), BP(J), C(J), D (J ) 3 E(J) , J = 1, . . . JM 

for all values of the spanwise index J. Allowance is made in the code for a 
term A T 3<f)/3t for which the value of the coefficient is stored in AP(J) and 
has been set to zero for all present uses. 

Initial and boundary conditions are required. To integrate this equation 
the initial condition used is usually the quasi-steady solution (i.e. , <j> xt = 0 
in eq. ( 1 ) ) . 

On the mean surface of the blade the flow tangency condition is expressed 
(cf. ref. 1 ) as: 


^ - (y + y cos t)f 1 (x) at z = 0 

At the innermost grid location, y m ± n > two boundary conditions can be 
used : 


1. A symmetry condition (equivalent to a flat tunnel wall in wing 

calculations) 


d± 

3y 


0 


(specified by setting JSYM = 1 in THREED) 
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2. A strip-theory condition (used for rotor blades or semi-infinite wings) 


3 2 <j> 

3y 2 


0 


In the far field a Dirichlet (<f> = 0) or a Newman (3<j>/3n = 0) condition has been 
used. The upstream boundary is usually taken as the uniform undisturbed flow 
(+ = 0). 


COORDINATE TRANSFORMATION AND THE CORRESPONDING MESH SYSTEM 


In order to treat a large class of planform shapes, a coordinate trans- 
formation is made prior to the discretization of the equation. This transfor- 
mation incorporates some one-dimensional stretching capabilities concentrating 
the mesh in regions of large gradients; in particular, near the surface of the 
blade, near the leading edge, and near the tip. The coordinate transformation 
is of the form, (x,y,z) -> 

where 


K - ?(x,y) 
n * n(y) 
e = c(z) 


Equation (1) now becomes: 


95 9 2 < j > / 95\ 2 3 2 4 > , P ,/ 3 C \ 3 9 / 3^\ 2 , L 95 9 ? . 

9x 9t 9? - B fc; VW 95 19?; + [ C 9x 9y + J 


9 2 <f> 
95 2 


9 2 (j> 
9C 2 


9 2 g 

9x 2 95 


, (r li M + 2 T> 21 la\ 

+ V° 9x 9y 2D 9y 9y,) 95 9n + ™9y/ 9t1 2 E \9z / 

o E , 95 9 2 5 / 9<J> \ 2 / 9 2 5 9 2 5\ 9<f> 9 2 n 94> 9 2 C 94> 

9x 9x 2 U/ + V C 9x9y + D 95 + D 9y 2 9 n + E 9z 2 9 C 


(2) 

The coefficients in equation (2) are the partial derivatives of the trans- 
formation. This form of the equation is called semi conservative; the metric 
coefficients are brought outside the 3/35, 3/3n» 3/3C symbols. It can be 
shown that, if the transformation is sufficiently regular, the jump conditions 
are preserved across a discontinuity. 

Computation is made of four first partial derivatives and five second 
partial derivatives. They are 3£/3x, 3£/3y, 3r|/3y, and 3c/3z (called in 
THREED XIX, XIY, YIY, and ZIZ, respectively) and 9 2 5/9x 2 , 9 2 5/9x9y, 9 2 5/9y 2 , 
9 2 Tl/9y 2 , and 9 2 ?/9z 2 (called in THREED XIX2 , XIXY, XIY2, YIY2, and ZIZ2, 
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respectively) . These quantities are computed at each interior mesh point by 
using finite difference approximations of the coefficients of the inverse 
transformation and the following identities: 


1 

1 

1 

0 


3x _3£ 
3? 3x 

3y 3p 
3q 3y 

3z 3g 
3C 3z 

3x 3g 

35 9y 


3x 3ri 
3ri 3y 


and , s imi lar ly , 

3 x 3 2 g 3 2 x /3|_\ 2 
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The derivatives 3x/3£, . . . are evaluated by finite differences at point 
i> j as 


/3x\ X i+i,.i X i-i,,j 

V951 , 2A? 

» J 


These expressions include second-order terms [0(A£ + An) 2 + 0(A£) 2 ]. The mesh 
is constructed in three steps. In the first step the locations of the span- 
wise stations are defined. The following stations are specified: 


y . 
mm 


innermost station on the blade, typically y . - 0.5 (y . is 

referred to as YN in THREED) rain min 

a special station on the blade (e.g., a kink in the planform) , 
typically y^ - 0.9 (y -> YC in THREED) 

the tip of the blade y _ = 1 (y -»• YD in THREED) 

u d 
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y outermost radial station, typically y - 1.5 ( y YX in 

max THREED) max X 

In addition to these real numbers the corresponding integers (JC < JD) must be 
defined. This determines how many stations for computation are located 
between 1 (y m i n ) and JC (y c ) , JC (y c ) and JD (y d ) , and JD (y d ) and JM (y max ) . 
The following analytical expressions are used to define the mesh stations: 

J > JD Y (J) = YX + cos^j y ~ ~~^ (YD - YX) 


J < JC Y (J) = YN + cos^|- — - — — ^(YC - YN) 

n - n 

JC < J < JD Y (J) = YC + — (YD - YC) 


where the variable r\ is defined between 0 and 1 by 


n 


J - 1 
JM - 1 


The planform of the blade then yields the locations x a and xf of the leading 
and trailing edges as functions of J. For this purpose, a piecewise analyti- 
cal representation is made of the planform. The chordwise coordinate trans- 
formation has no radial dependence for all points beyond the tip. In THREED 
x a and Xf are called XA(J) and XF(J). 

The second step in mesh construction, in the chordwise direction, is 
defining 

Xmin upstream boundary, typically - -8 (x^ -> XN) 

^ax downstream boundary, typically x max - 6 (x max -*■ XX) 


and the indices IA < IF which determine how many stations are located 
between 1 (x^n) and IA (x a ) , IA (x a ) and IF (x f ) , and IF (xf) and IM (x^x) . 
Similar analytical expressions are used to define the mesh stations in x: 

I > IF X(I,J) = XX + cosfe x _ A xF(J) - XX] 


I < IA X(I,J) = XN + cosfe j[XA(J) - XN] 


IA < I < IF 


X(I,J) = XA(J) + 



[XF (J) - XA(J) ] 
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where the variable E. is defined between -1 and 1 by 


5 = -1 


2(1 - 1 ) 

IM - 1 


In the third step in the vertical direction the following are defined: 


z min lower boundary, typically z min - -3 (z m i n -> ZN) 

z max upper boundary, typically z may - 3 (z max -v ZX) 


and the indices KU = KO + 1 which determine how many stations are located 
between 1 (z m j_ n ) and KO (nearest to the lower surface of the blade) , and KU 
(nearest to the upper surface of the blade) and KM (z max ) . The mesh stations 
in z are defined by using the analytical expressions: 


K > KO Z(K) = ZX - cos(|- ?)ZX 

K < KO Z(K) = ZN - cos(| ^ZN 


where £ 


is defined between -1 

C 


and 1 by 


= -1 + 


2(k - 1) 

KM - 1 


The mesh dimensions in the code are set up to allow for maximums of IM = 64, 
JM = 32, and KM = 32. 


FINITE DIFFERENCE SCHEME 


In equation (1), the nonlinear term (d /dx) [B (d<p/3x) + B T ( d<j>/dx) 2 ] , which 
is often written nonconservatively as V(J> XX , is responsible for the mixed 
character of the flow. It is well established that a mixed scheme must be 
used for the nonlinear flux discretization (refs. 2 and 3), given as follows 
for a uniform mesh spacing: 


Let 


V. = B + 2B* 
x 


i+i 


- *±-i 


2Ax 


In the following four cases to be considered the nonlinear term is discretized 
e.g.. 

Case 1 V_^ > 0 - 0 (subsonic point) 

*i+i - 2 *i + *±-i 

Discretization: V, t : — 

i Ax 2 

\ 

(the indices j and k, which are invariant, are not indicated). 


7 



Case 2 < 0 V^_ x < ® (Supersonic point) 


Discretization: 



24> i-l + *1-2 
Ax 2 


Case 3 V. < 0 V. > 0 
1 1-1 


(sonic point) 


Discretization: 



2<j> 


1-1 


+ * 


i-2 


Ax" 


Case 4 V. > 0 V. <0 (shock point) 
i 1-1 


Discretization: 


<f>. , - 2* . + <f>. 

V 1 1 — - — - — + v 

i Ax 2 i- 1 


(j) . — 2cf) , + <b . 

l l-i 1-2 

Ax 2 


In contrast to most small disturbance codes (typified by refs. 4 and 5), 
the discretization of the sonic point (case 3) eliminates some spurious oscil- 
lations that appear when the sonic line is located close to the leading edge 
of a blunt airfoil in a region where the flow experiences a rapid acceleration. 
It can be shown that the discretization that is proposed here is consistent 
with the equation, but it is not strictly conservative. However, the error of 
conservation is small, and not larger than 0(Ax). The shock-point discre- 
tization, however, ensures conservation of mass at the shock point. 


The next term in equation (1) is the cross-derivative term. This term is 
small inboard where the flow is subsonic and two-dimensional. However, for 
large advance ratios (y - 0.5) and for values of azimuth and radius where the 
transonic flow has a large radial component, its effects cannot be neglected. 
In fact, in these cases the cross-derivative term, which is usually treated 
explicitly (i.e., always at the previous time level (ref. 5)), has a destabi- 
lizing effect and can strongly reduce the time step required for maintaining 
overall stability. 


For values of C > 0, corresponding to a negative sweep angle, the cross- 
derivative term is discretized as: 


c . 
3 


i»j+i 


|> . . - <|> . . , + <p . 

1,3 Y i-i,J+l i-i >3 

Ax Ay 


For values of C < 0, corresponding to a positive sweep angle, the following 
discretization is used: 


c . 
3 


The schemes that are presented for uniform mesh spacing extend readily 
to the mesh obtained from the coordinate transformation. The coefficient of 
the cross-derivative is now 
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c 


m . 2t> la 

3x 3y 3y 3y 

For discretization of the term c (3£/3x) (3£/3y) (3 2 cf>/3£ 2 ) a centered scheme is 
used at all points: 

3g. . 3£. . (f) . - 2(J>. + (f) . 

i»J Y i+i T i 1-1 

j 3x 3y A £ 2 

SOLUTION ALGORITHM 


The time-accurate integration is obtained by using an Alternate Direction 
Implicit (ADI) scheme that breaks the three-dimensional problem into three one- 
dimensional problems in each coordinate direction. The advantage of this 
scheme is its inherent stability, at least in the case of a linear equation, 
regardless of the local Courant number. Indeed, when solving a complicated 
problem, in a mesh where the cell sizes may vary by one or more orders of mag- 
nitudes, it would be very time-consuming to limit the time step to satisfy the 
Courant number associated with the smallest cell. 

However, since the equation being solved is nonlinear, there is a prac- 
tical limitation which can be associated with vortex shedding in lifting 
cases or configurations with shock motion. This means that the cell sizes 
must never be so small on the airfoil surface and near the trailing edge 
that the allowable time step for maintaining stability is unnecessarily 
limited. 


A Crank-Nicholson averaging between the levels n and n + 1 is used 
since it can be shown on the linearized equation that a stable scheme results. 

The three steps of the ADI-Crank-Nicholson algorithm are as follows: 

Step 1 


15 A. (± - A m 1 /HN VaV , 3 2 A . R ,/15\ 
9x 95 \ At / 2 \3x/ \3 £ 2 9 £ 2 / \9x/ 

+ ( C If 3y 

+ D (Al\ 9V 1 . + e /9JL\ 

\9y/ 3 n 2 \9z/ 


3 

95 


3<j> n 3j>_ 

95 9? 


D /3£\ /9fj n . 9V 

2 \9y/ \ 35 2 3£ 2 / 


+ 2D — 


la \ 9 2 <f> n 

3y 3y/ 3£3 ti 


+ c 


3g 35. 3 2 <}> n 
3x 9y 9^2 


2 v 


1 B 3 2 5 
2 3x 2 


\ 95 95/ 


9r 

2 i 

9x 3x 2 


+ 2 B' 15 1VM 


n ~ 


95 95 


M. + D 3 g 


9y 


/V + ii\ 

\ 95 95/ 


+ D 


3 2 n 3(j> n 
3y 2 9n 


n 


+ E 


3 2 g 3£ 
3 z 2 3? 


+ c 


_3£ 


n 


9 2 5 

3x3y 3g 
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It should be noted that the underlined terms are treated explicitly. However, 
an implicit scheme can easily be devised based on switching from a centered 
difference approximation to 9 2 cf>/3£ 2 when (c3£/3x) (3£/3y) > 0, and to an 
upwind difference approximation when (c3£/3x) (3£/3y) < 0. Test results for a 
swept tip showed very little difference between the implicit treatment and the 
explicit scheme given previously. 

Step 2 

A H _L ( l - A = I ( e li 9n + 2n K 3n\ / 9 2 A _ 9 2 <j> n \ 

9x 95 \ At / 2 V 9y 9y 9y/\ 959n 959n / 

D /in\ 2 /9fi 9 2 A, D 3 2 n /s* Scf> n \ 

+ 2 \9y/ (a,-, 2 ^ 2 / 2 3 y 2 y9n 9n / 

Step 3 

.95 _9_ /<j) n+1 - = E / 9g\ 2 / 9 2 j» n+1 _ aVV E 9f£ /9j» n+1 _ 9(j> n \ 

9x 95 \ At / 2 V9z/ \ 9C 2 9 ? 2 / 2 8z 2\ 3 C 3 ? / 

After these equations are discretized according to the method discussed in the 
Finite Difference Scheme section, the algebraic system is inverted by using a 
tridiagonal or quadradiagonal direct solver. Particular attention is given, 
when defining the finite difference analogues, to ensure that the main diagonal 
term could be chosen as pivot in the elimination process. All the terms, which 
are treated implicitly, contribute to the main diagonal with the same sign. 

Each complete time step requires approximately 2.5 sec of CPU time of the 
CDC 7600 computer. A rectangular blade computation requires approximately 
half an hour of total run time. For swept tips, where there is a more severe 
time-step limitation, the total run time is an hour. The corresponding mesh 
is composed of approximately 35,000 nodes. 


RESULTS 


Some three-dimensional steady (hover) flows are simulated for three blade 
geometries: a rectangular blade, a swept-tip blade, and a blade combining a 

swept and parabolic tip (fig, 1). The pressure distributions are presented 
for three sections of the blades in figures 2(a-c) for blade A, figures 2(d-f) 
for blade B, and figures 2(g-i) for blade C. As can be seen, the effect of 
sweep is favorable inboard. The shock waves either are weakened or disappear 
on blade C. Close to the tip, however, the opposite trend seems to occur, with 
blade C experiencing the largest supersonic pocket. The global effect is 
favorable for the swept tips in hover. 

Three-dimensional unsteady flows past a rectangular blade of aspect ratio 
/R = 7, have been computed at Mach numbers M = 0.6 and advance ratios (y) of 
0.45, 0.5, and 0.55. The blade has no twist and is equipped with a symmetric 
NACA 00XX profile of varying thickness along the span. ONERA experimental 
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data for the same rotors and test conditions are available for comparison 
(ref. 4). Figure 3 shows the radial stations for the experimental pressure 
measurements. The corresponding results are shown in figures 4(a-c) for the 
azimuth of 60° and in figures 4(d-f) for the azimuth of 120° at the lowest 
advance ratio. For the advance ratio of y = 0.5, the results are presented 
in figures 4(g-l). Figures 4(m-r) show results for y = 0.55 for the same 
two azimuth angles of 60° and 120°. Also plotted in these figures are the 
quasi-steady results, which correspond to 3 2 (|)/3t3x = 0. As can be seen, the 
unsteady results agree better with the experimental results, indicating a non- 
negligible unsteady term 8 2 (f>/3t3x. Furthermore, a comparison of the quasi- 
steady solutions at azimuth angles of 60° and 120° exhibits the influence of 
the cross-derivative 3 2 <f)/3x3y, which increases toward the tip. 


THREED CODE 


THREED has been coded in FORTRAN by using only standard statements. In 
its present form it is adapted to the CDC 7600 computer of the Ames Research 
Center, NASA, Moffett Field, California. The Small Core Memory (SCM) length 
is 27,257 decimal words and the Large Core Memory (LCM) length is 131,072 deci- 
mal words. THREED is divided into one main program and four subroutines: 

SUBROUTINE MESH defines the mesh and computes the metric coefficients 
SUBROUTINE SLOPE computes the slope of the blade at each point 
SUBROUTINE POT integrates the potential equation 

SUBROUTINE CP computes the pressure coefficient on the blade 

The data as they are read in and printed out are shown in appendix A. The 
values shown correspond to the results plotted in figure 4(g-l). A short 
explanation of the parameters as well as the notation in THREED follows: 

i) HM = 0.6 Mach number 

ALPHO = 0 mean incidence of the sinusoidal motion, deg 

DALPH = 0 amplitude of the incidence variation, deg 

IR0T = 1 rotating blade case 

IR0T = 0 fixed blade or wing case 

AV =0.5 advance ratio 

GM = 1.4 ratio of specific heats 

ii) NSTP = 601 number of time steps 

ITER = 400 number of iteration steps 

DTN = 0.0001 minimum time-step size in the relaxation process (rad) 

DTX = 0.01 maximum time-step size in the relaxation process (rad) 

NM0D = 8 number of elements in the time-step sequence based on 

DTN and DTX 

TI = -1.5708 initial time (rad) 

NPR = 100 time step at which results are printed 

iii) YC = 0.9 location of the kink or a special span location 

YD = 1 tip of the blade 
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iv) 


v) 


vi) 


KSG = 1 

KSG = 0 
DEL = 0.12 
AR = 7 

JYSM = 0 

JYSM = 1 
KSYM = 1 
KSYM = 0 
KGRAD = 1 
KGRAD = 0 

IM = 64 
JM = 32 
KM = 32 

IA = 18 
IF = 48 
JC = 11 
JD = 21 
KO = 16 
KU = 17 


indicates that the blade is symmetric with respect to 
z = 0 

indicates that the blade is not symmetric 

thickness of the basic profile as defined subsequently 

aspect ratio of the blade 

indicates that a strip-theory condition is used at the 
root 

indicates that a symmetry condition is used at the root 
there is a lower-upper symmetry 
there is no lower-upper symmetry 
a Neuman boundary condition is used 
a Dirichlet boundary condition is used 

number of mesh points in the E, direction ( <64 ) 

number of mesh points in the n direction ( <32 ) 

number of mesh points in the £ direction ( <32 ) 

leading-edge index (IA < IF) 
trailing-edge index (IF < 64) 
kink station index (JC < JD) 
tip station (JD < 32) 
lower-surface index 

upper-surface index (KU = KO + 1 < 32) 


vii) 


XN = -8 location of most upstream surface 

XX = 6 location of most downstream surface 

YN = 0.5 location of innermost surface 

YX = 1.5849 location of outermost surface 
ZN = -3 location of most bottom surface 

ZX = 3 location of most tip surface 


viii) The basic airfoil is defined by two sets of values representing the 
abscissas and the ordinates of points on the profile. The maximum number of 
points is 101. The points are distributed in sequential order around the air- 
foil, starting from the trailing edge describing the upper surface, then the 
lower surface, and then back to the trailing edge. 


NP = 101 total number of points (<101) 

(NP -1- l)/2 must correspond to the leading edge 


When the geometry changes rapidly it is preferable to concentrate the 
points near the leading and trailing edges in order to ensure the best pos- 
sible accuracy for linear interpolation. The coordinate profile at the mesh 
location N is 


XP (N) 
ZP (N) 


abscissa 

ordinate 


of point N 


The planform geometry of the blade is defined by piecewise-analytic formulas 
in the subroutine MESH. An example is given in appendix B for a swept tip. 
The functions XA(J) and XF(J) are defined in the loop starting with 
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77 


DO 6 J = 1, JM 


and ending with 

87 6 CONTINUE 

as shown in the box. 


CONCLUSIONS 


A finite difference code for predicting the high-speed flow over an 
advancing helicopter rotor is presented. The code solves the low-frequency 
transonic small disturbance equation and is suitable for modeling the effects 
of three-dimensional advancing blade unsteadiness. This work was inspired by 
a similar method developed by F. X. Caradonna (ref. 5). However, the computer 
code THREED incorporates some important new features, especially the capabil- 
ity for treating nonrectangular blade tips. Computed results show good 
agreement with experimental blade pressure data and illustrate some of the 
effects of varying the rotor planform. The flow unsteadiness is shown to be 
an indispensible part of a transonic solution. It is also shown that close to 
the tip at high advance ratio, cross-flow effects can significantly affect the 
solution. 


Ames Research Center 

National Aeronautics and Space Administration 
and 

Aeromechanics Laboratory 

AVRADCOM Research and Technology Laboratories 

Moffett Field, Calif. 94035, April 10, 1980 
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APPENDIX A 



SAMPLE OF 

DATA AS READ IN AND 

PRINTED OUT 

OF CDC 7600 

COMPUTER 


MACH 

NUMBER 

.6 

MEAN 

INCIDENCE 

0. 

INCIDENCE 

VARIATION 

0. 

ROTATION 

Y/N 

1 

ADVANCE 

RATIO 

.5 

HFAT 

RATIO 

1.4 



NO. TIME 
STEPS 
601 

NO. STEPS 
RELAXATION 
40 

MIN STEP 
RELAX 
.0001 

MAX STEP 
RELAX 
.01 

RELAXATION 

CYCLE 

8 

INITIAL 

TIME 

-1.5708 

FINAL 

TIME 

1.5708 

IMPR. 

STEPS 

100 

SPECIAL 

SPAN 

LOCATION 

.9 

TIP 

LOCATION 

1. 

GEOM. 

SYMMETRY 

Y/N 

1 

BASIC 

THICKNESS 

.12 

ASPECT 

RATIO 

7. 




ROOT SYM 
CONDI TION 
Y/N 
0 

UP-LO 

SYMMETRY 

Y/N 

1 

LATERAL 

GRADIENT 

Y/N 

1 







NO • X 

NO. Y 

NO • Z 

MESH 

MESH 

MESH 

64 

32 

32 


LEADING 

TRAILING 

SPECIAL 

TIP 

LOWER 

UPPER 

EDGE 

EDGE 

SPAN NO. 

NO. 

SURFACE NO. 

S. NO. 

18 

48 

11 

21 

16 

17 


MIN.X 

MAX.X 

MIN. Y 

MAX .Y 

MIN.Z 

MAX.Z 

SURFACE 

SURFACE 

SURFACE 

SURFACE 

SURFACE 

SURFACE 

-8.0 

6.0 

.5 

1.5849 

-3. ■ 

3. 
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BASIC PROFILE 
NO. POINTS INDEX 
101 

1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


ABSC I SSA 

ORDINATE 

1.0000 

.0013 

.9990 

.0014 

.9961 

.0018 

.9911 

.0025 

.9843 

.0034 

.9755 

. 0046 

.9649 

.0061 

.9524 

.0077 

.9382 

.0096 

.9222 

.0117 

.9045 

.0139 

.8853 

.0163 

. 8645 

.0188 

.8423 

.0214 

.8187 

.0241 

.7939 

.0269 

. 7679 

.0297 

.7409 

.0325 

.7129 

.0354 

.6841 

.0382 

.6545 

.0409 

.6243 

.0436 

.5937 

.0461 

. 5627 

.0486 

.5314 

. 0500 

.5000 

.0529 

.4686 

.0548 

.4373 

.0564 

.4063 

.0578 

.3757 

.0588 

.3455 

.0596 

.3159 

.0600 

.2871 

.0600 

.2591 

.0596 

.2321 

.0589 

.2061 

.0577 

.1813 

.0562 

.1577 

.0542 

.1355 

.0519 

.1147 

.0491 

.0955 

.0460 

.0778 

.0426 



BASIC PROFILE 

NO. POINTS INDEX ABSCISSA ORDINATE 
101 


43 

.0618 

.0389 

44 

.0476 

.0348 

45 

.0351 

.0305 

46 

.0245 

.0259 

47 

.0157 

.0211 

48 

.0089 

.0161 

49 

.0039 

.0109 

50 

.0010 

.0055 

51 

.0000 

.0000 

52 

.0010 

-. 0055 

53 

.0039 

-.0109 

64 

.0089 

-.0161 

55 

.0157 

-.0211 

56 

.0245 

-.0259 

57 

.0351 

- .0305 

58 

.0476 

-. 0348 

59 

.0618 

-.0389 

60 

.0778 

-.0426 

61 

. 0 Q 5 5 

- . 0460 

62 

.1147 

-.0491 

63 

.1355 

-.0519 

64 

.1677 

-.0542 

65 

.1813 

-.0562 

66 

.2061 

-.0577 

67 

.2321 

- .0589 

68 

.2591 

-.0596 

69 

.2871 

-.0600 

70 

.3159 

-.0600 

71 

. 3455 

-.0596 

72 

.3757 

-.0588 

73 

.4063 

-.0578 

74 

.4373 

— .05 64 

75 

.4686 

-.0548 

76 

. 5000 

-.0529 

77 

. 5314 

-.0509 

78 

. 5627 

-. 0486 

79 

. 5937 

-.0461 

80 

.6243 

-.0436 

81 

.6545 

-.040° 

82 

.6841 

-.0382 

83 

. 7129 

-.0354 

84 

.7409 

-. 0325 

85 

.7679 

-.0297 

86 

.7939 

-. 0269 

87 

.8187 

-.0241 

88 

.8423 

-.0214 


16 



BASIC PROFILE 
NO. POINTS INDEX 
101 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 


ABSCISSA ORDINATE 


.8645 

-.0188 

.8853 

-.0163 

.9045 

-.0139 

.9222 

-.0117 

.9382 

-.0096 

.9524 

-.0077 

.9649 

-.0061 

.9755 

-.0046 

.9843 

-.0034 

.9911 

-.0025 

.9961 

-.0018 

.9990 

-.0014 

1.0000 

-.0013 


APPENDIX B 


SUBROUTINE MESH 


DUM= ( Y T J-YIC) / ( Y 10 — VIC) 

Y ( J )=YC+01JM* ( Yf)-YC ) 

c;n Tn a 

2 CONTINUE 

OI.H = - . 5* P I * ( Y I J- Y I C ) / Y I C 
DllM=cnS ( HUM ) 

Y ( J )=YN+OUM* ( YC-YN ) 

on Tn a 

3 CONTINUE 

DtJM= . 5*PT* ( Y I J-YI 0) / ( 1 .-Ylf) ) 
DiiYi=cns ( hum) 

Y { J ) =YX+niJM* ( YO-YX ) 

4 CONTINUE 

I F ( JSYM.ME. 1 ) Y ( 1 ) = Y ( 2 ) 

WR I TF ( 6, 1000 ) 
nn 5 J = 1 t JM 

WRITE (6» 1001 ) JrYI (J) »Y(J) 

5 CnNTINUE 

c *****p LAM FORM EQUATION XA ( J ) * XF ( J ) 
on 6 J=1»JM 

Y J = Y ( J ) 

XA ( J )=0. 

XF( J)=l. 

IF (J.LF.JC) GU TO 6 
XA( J ) = 3. 531767*( YJ-YC) 

XF( J)=1.+XA( J) 

IF(J.LE.JO) GU TG 6 
XA( J ) = X A ( J-l ) 

XF< J ) = XF ( J-l ) 

6 CONTINUE _____ 

X A ( 1 ) = X A ( 2 ) 

XF ( 1 )=XF (2) 

WR I TF (6, 1002 ) 
on 7 J=1*JN 

WRITE (6» 1001) J* XA( J > »XE { J) 

7 CONTINUE 

XI I=-1.-DXI 
on 8 I = 1 , I M 
xi i=xi i+nxi 

XI (I ) = X I I 

8 CONTINUE 
1=0 
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M = 0 



Figure 1.- Simulations of three-dimensional steady (hover) flows for: A - a 

rectangular blade; B - a swept-tip blade; and C - a combination of swept- 
and parabolic -tip blade. 
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JU = 0.45 , 0.50, and 0.55 


/OR = 7 



RECTANGULAR BLADE 



Figure 3.- Three-dimensional unsteady problem for forward flight. 
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Figure 4.- Continued 
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Figure 4.- Concluded 
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